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Abstract 

In a companion report, we have derived a method for finding the impedance at high frequencies 
of vacuum chamber transitions that are short compared to the catch-up distance, in a frequency 
regime that — in analogy to geometric optics for light — we call the optical regime. In this report 
we apply the method to various non-axisymmetric geometries such as irises/short collimators in a 
beam pipe, step-in transitions, step-out transitions, and more complicated transitions of practical 
importance. Most of our results are analytical, with a few given in terms of a simple one dimensional 
integral. Our results are compared to wakefield simulations with the time-domain, finite-difference 
program ECHO, and excellent agreement is found. 
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I. INTRODUCTION 



In many current and future accelerator projects short, intense bunches of charged particles 
are transported through vacuum chambers that include objects such as transitions, irises, 
and collimators. For example, in the beam delivery system of the International Linear 
Collider (ILC) 3 nC, 300 /im-long bunches encounter many collimators on their way to 
the interaction region [1]. Or, in the undulator region of the Linac Coherent Light Source 
(LCLS) a 1 nC, 20 /im-long bunch passes by square-to-round transitions, bpm's, and other 
changes in chamber geometry [5] . Wakefields generated by such changes in vacuum chamber 
geometry can negatively affect the beam emittance and ultimately the performance of an 
accelerator. Numerically obtaining the strength of the wakefields or, equivalently, of the 
impedances for short bunches in such vacuum chamber objects can be difficult, in particular 
when the object is long and non-cylindrically symmetric. 

In a recent paper [3] we developed a method to solve such seemingly difficult problems in 
a so-called optical approximation. This method is valid in the limit of high frequencies and 
reduces the calculation of the impedance to the two dimensional integration of potential 
functions. In this report we make use of this method to work out solutions to a selection 
of 3D (shorthand for non-cylindrically symmetric) geometries that can be encountered in 
today's accelerators. 

The geometry of the problems to be considered is, in general, of the type sketched in 
Fig. [I] An in-going beam pipe (region A with cross-section profile Sa) is followed by a 
short transition (the gap or aperture region G with cross-section Sg) and ends in an out- 
going pipe (region B with cross-section Sb)- We limit consideration to cases where Sg is 
contained within the intersection of Sa and Sb [Sg ^= (Sa H Sb)]- The transition need 
not be smooth like the one in the figure. A speed of light beam passes through the three 
regions on a straight line that we call the design orbit. At first there is no assumption on 
transverse symmetry in the three region. However, we do assume that the axes of regions 
A and B are parallel to the design orbit. Note that if Sg = Sb the structure is called a 
step-in transition (by which we mean a short, inward transition), if Sg = Sa it is a step-out 
transition (a short, outward transition). If Sa = Sb, with the aperture of Sg smaller than 
the beam pipes, then the structure is an iris (a metallic diaphragm with a hole or slot) or 
short collimator in a beam pipe. 




FIG. 1: Sketch of a generalized 3D transition, showing the regions A, G, and B, with cross-sections 
Sa, Sg, and Sb, respectively. 

We are interested in finding the high frequency longitudinal or transverse impedance of a 
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general transition such as is sketched in Fig. [T] The impedance regime can be described as 
the regime of geometric optics. This regime is applicable provided that (1) the frequency is 
high and (2) the transition is short. By condition (1), we mean that the frequency uj 3> c/g, 
with c the speed of light and g the minimum aperture of the structure. If the transition is 
tapered, we require further that to S> c/(g9), with 9 the taper angle. By condition (2), we 
mean that the length of the transition is short compared to the catch- up distance, I = g 2 u/c. 
Note that in the optical regime the high frequency longitudinal impedance of a structure 
is resistive, with the impedance Z\\ real and independent of frequency. The transverse 
impedance Z± is also real and depends on frequency as u . 

The method of Ref. [3] can be used to find the impedance, for example, for an iris in a 
beam pipe, a step-in or step-out transition, and a short collimator. Also, there are objects 
that we call long collimators; i.e. collimators that have a fixed (minimum) aperture over a 
length that is long compared to the catch-up distance £, with transitions at the front and 
back ends. For long collimators our methods also apply, provided that the two transitions 
satisfy the above two conditions. For such a collimator the impedance is the sum of the 
impedances of its two transitions. Note, however, that for intermediate length collimators, 
those of length comparable to the catch-up distance, our methods do not apply. 

In Ref. [3] the impedance of a general transition in the optical regime is developed follow- 
ing a systematic approach. Earlier work on the subject was focused on specific geometries 
and often was rather informal in nature. Balakin and Novokhatski were among the first to 
address the question of impedance in the optical regime |4J. Hcifcts and Khcifcts studied 
the longitudinal impedance of round step-in and step-out transitions (in the optical regime) 
numerically by field matching [5J. For round, long collimators the dipole mode impedance 
was obtained rigorously in Ref. [B], and the higher azimuthal mode impedances, more infor- 
mally, in Ref. [7]- The impedance of short and long, round and 3D collimators was studied 
numerically in Refs. [8], it was found, for example, that the impedances of short and 
long collimators, in general, differ significantly in amplitude (by a factor ~ 2 in the round, 
transverse case). 

The impedances that we obtain in this report following our method will also be compared 
to numerical results obtained by the computer program ECHO [10] . This program solves 
Maxwell's equations to find the wakefield of an ultra-relativistic Gaussian bunch within 
(perfectly conducting) metallic boundaries of fully 3D geometry. From the wakefields of a 
sufficiently short bunch the impedance in the optical regime can be obtained directly. In 
the present report discussions of ECHO calculations will be brief; their main purpose is to 
confirm our results and to give us confidence in our method. 

This report is organized as follows: In Section II the method of calculation, derived in 
the companion paper, Ref. [3], is presented. At the end of this section some details of the 
ECHO simulations will be given. The heart of the present report, however, is the next 
four sections where our method is applied to 3D transition examples. The problems that 
we solve are example irises or short collimators within a beam pipe (Section III), step-in 
transitions (Section IV), step-out transitions (Section V), and more complicated transitions 
(Section VI). In these sections our results will be compared briefly to ECHO simulation 
results. Section VII gives the conclusions. In the Appendix specific limits for the impedance 
of an elliptical step-out transition are derived. Note that although we follow the method 
of Ref. [3] in this report, the notation used here is not exactly the same. Note also that 
Gaussian units will be used throughout; to convert impedances to MKS units, one multiplies 
by the factor Zqc/Att, with Zq = 377 fl. 
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II. IMPEDANCE CALCULATIONS 



Consider a general (transversely non-symmetric) transition. Let the design orbit follow 
the z-axis, with the particle motion in the +z direction. The transverse impedance of a 
general transition consists of monopole, dipole, and quadrupole components (with respect 
to the design orbit) that involve tensors, and the method of Ref. [3] can, in principle, deal 
with such problems. However, to simplify, in this report we will limit consideration to 
geometries for which the design orbit lies on a vertical symmetry plane of the boundaries, 
defined by (horizontal coordinate) x — 0. Then the transverse (vertical, y) impedance (with 
respect to the reference trajectory) can be written as 

Zj_,tot = Z±^ m + yiZ± t d + V2Zx,q , (1) 

where the three terms are called the (transverse) monopole, dipole, and quadrupole contri- 
butions; where y\ is a small offset of particle 1 (the leading particle), and t/2 a small offset 
of particle 2 (the trailing particle, see Fig. |2j. Note that the horizontal impedance has an 
equation equivalent to Eq. (flj), with the quadrupole term equal to (— Z± i9 ). 




FIG. 2: Sketch of a transition showing the relationship of the design trajectory and the vertical 
offsets of particle 1, j/i, and particle 2, j/2 (the particles are indicated by red dots). For the 
calculations the offsets yi, j/2, are assumed to be small. 

For most of the examples of this paper, the system has also a horizontal symmetry plane 
and the design orbit lies in this plane. In such a case Z±_ m = 0; if, in addition, y\ = y% = yo 
we can define a normalized total impedance 

Z± = ±,fot = Z± t d + Z± <q . (2) 

By convention, this is the normal definition of total transverse impedance for bi-symmetric 
problems. (In this paper, however, we will calculate the individual terms and Zj_ iQ as 
well as Z± for bi-symmetric problems.) The transverse wake at position s within a bunch 
(s < is toward the head), normalized to the offset yo, in the optical regime is given by 

W±(s) = ± ^ Z ± (tj)\ z (cj)e-^ s / c dw = (ujZ ± ) ^ X z (s') ds' , (3) 

with \ z (uS) the Fourier transform of the line charge density A z (s). To arrive at the last 
expression in this equation we have used the fact that, in the optical regime, Z^_ oc w . 
Thus the quantity (loZ±) is independent of frequency. The result is proportional to the 
integral of the line charge density (and in the longitudinal case, it is proportional to the 
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line charge density itself). The kick factor n±, the average kick experienced by the beam 
per unit charge per unit offset, is the integral of Wj_(s) when weighted by the longitudinal 
charge density. In the optical regime the kick factor is simply given by the constant 



Kj_ 



(4) 



In Ref. [3] using the Panofsky-Wenzel theorem [TT], the transverse monopole, dipole, and 
quadrupole terms are related to longitudinal monopole, dipole, and quadrupole impedances: 



Z±,m — 



Z-L,d = —Z\\,d 



Z^.q — Z\\ 



(5) 



The longitudinal impedances can, in turn, be obtained from integrals involving the Green 
functions to Poisson's equation (the potentials) in regions A and B of the transition of 
interest. In the case of the monopole part 
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V<^ m ,B • Vv5d,s dS 
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^7<Pm,A ■ Vyrf.fl dS 



s G 



(6) 



with the integrals taken over the cross-section areas Sb and Sq, respectively. Here the 
monopole and quad potentials in region A are given by the solutions to 



-4irS(y)6{x) , V 2 ^ = 4n5'(y)6(x) 



(7) 



with boundary conditions ip m ,A = 0, tfd,A — 0, on metallic boundary C'a that encloses Sa- 
Similar equations hold for ip m ,B and ipd,B of region B. 
In the case of the dipole part of the impedance 



Z\ 



,d 
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For the quad part of impedance 
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/ ^(Pm,B ■ ^<fiq,B dS - V(p m ,A ■ V<fiq, 
J Sb J Sg 



■dS 



The quadrupole potential in region B is given by the solution to 

VVb = -2*6"{y)&{x) , 



(8) 



(9) 



(10) 



with ipq,B = on boundary Cb- Finally, the longitudinal impedance on the reference 
trajectory is given by: 



,lo7ig 
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2ttc 



/ (^fm,B) 2 dS- Vip m ,A ■ V<£ m , 

J Sb J Sg 



■dS 



(11) 



To obtain the needed potentials ip m , ipd, ip q , we begin with the Green function solution 



to 



W 2 G(x, y, y ) = -4nS{x)S{y - y ) , 



(12) 



where G = on the boundary C of the surface of interest S. Here the source particle is 
located at x = 0, y — yo- For region S, tp m (x, y) — G[x, y, 0), 



<Pd(x,y) 



d 

dy 



G(x,y,y ) 



Vo=0 



dy; 



G(x,y,y ) 



(13) 



y =o 
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As was shown in Ref. [3J, with the help of Green's first identity fT2"] . 

/ (4>\7 2 ip + V0 • dS = I (/m-V^dl , (14) 

the surface integrals involved in the impedance equations can be converted to line integrals. 
In this formula <j> and ip are functions that can be differentiated twice; S is a surface over 
which a surface integration is performed; C is the contour that encloses S and over which a 
line integral is performed; and n is a unit normal vector pointing outward from the contour. 
This device is used throughout the examples of this report. 



Numerical Comparisons 

A great number of ECHO simulations were performed to test and verify our results. ECHO 
is a 3D, time-domain finite difference program that calculates wakefields generated by an 
ultra-relativistic bunch passing through the structure. ECHO has two features that make 
these 3D, optical regime simulations tractable. (1) A method to reduce the so-called "mesh 
dispersion" — errors generated in time-domain mesh programs that are especially difficult 
to deal with for the combination of short bunches and long structures. (2) An indirect 
method of calculating wakes in 3D structures that eliminates long downstream beam pipes, 
which would be needed in the case of direct wake calculation with very short bunches. Since 
the ECHO simulations play a supporting role in the present report, other than comparing 
results, we will only discuss them briefly. For a more detailed report on ECHO, see e.g. 
Ref. [ID]. 

A few comments on the parameters used in the ECHO simulations: The bunches in the 
simulations are Gaussian. Typically we choose a z — g/10, with g the minimum vertical half 
aperture of the iris or transition. For a flat iris or beam pipe we take the structure half-width 
to = lOg. (In this report we use the word flat to mean "having a rectangular cross-section 
with a small height to width ratio." ) For a small iris in a beam pipe or a transition opening 
into a large beam pipe, the beam pipe is square with half-height lOg. We take mesh sizes to 
be |c 2 , §Cz, in the longitudinal, horizontal, vertical directions, respectively. Irises have 
the thickness of one longitudinal mesh size. For bi-symmetric structures ECHO calculates 
the dipole and quadrupole components of the wakes independently. 

To demonstrate the validity of the optical approximation, and how one moves out of the 
optical regime as the bunch length increases, we consider the problem of a thin round iris 
of radius g in a round beam pipe of radius b. We choose b/g = 4. We perform ECHO 
calculations for bunch lengths in the range <r z fg = 0.02 to 4. In Fig. [3] we plot twice the 
kick factor of our results as function of bunch length. According to Eq. |4j in the optical 
regime 2k±^ — ojZ±_ , which is a number independent of bunch length or frequency We see 
that the result of our numerical calculation is fairly constant up until <J z /g ~ 0.2, after 
which the validity of the optical approximation starts to break down. In more detail, we 
give in Fig. [4] the actual wake functions for four example bunch lengths as obtained by 
ECHO. The (Gaussian) bunch shape, with the head to the left, is indicated by the black 
dashes. The analytical wake of a Gaussian bunch in a round iris in the optical regime, 
W± = [1 +erf(s/V2a z ))/g 2 , with erf(x) the error function, is also given (the red dashes). 
We see that, although the weighted average of the wake agrees well for a z /g ~ 0.2, good 
agreement between the wakes over ±4cr z is not obtained until a z j g ~ 0.04. 
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FIG. 3: Twice the kick factor k± of a thin round iris in a beam pipe vs. rms bunch length cr z , as 
obtained by ECHO. The iris radius is g. 
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FIG. 4: Transverse wake of a thin round iris in a beam pipe as obtained by ECHO. Results for bunch 
lengths over iris radius cr z /g — 0.04, 0.16, 0.4, and 0.8 are given. The (Gaussian) bunch shape, with 
the head to the left, is indicated by the black dashes. The analytical wake of a Gaussian bunch in 
a round iris in the optical regime, W± = [1 + eri(s/\/2a z )]/g 2 , with erf(x) the error function, is 
also given (the red dashes). 



III. IRIS/SHORT COLLIMATOR IN BEAM PIPE 

We begin by performing calculations of examples with thin irises or short collimators in a 
beam pipe (for an example, see Fig. ISJ. In these cases the solution is given by a particularly 
simple "clipping" type calculation of the energy impinging on the iris wall. 

For these examples the beam pipes in Region A and B are identical, i.e. Sb = Sa- Thus 
Eqs. (§, (§, (§, give 

Z±,d= ^- / (Vvd, B ) 2 dS , (15) 
Z±, q = — / V</? m .B • Vp q ,B dS . (16) 

™ JSb-Sg 



7 



FIG. 5: An iris in a beam pipe. 



Both equations involve integrals over the metallic surface of the iris. Using Green's first 
identity and considering the fact that Vy 2 = in any region that precludes the axis, we 



can write Eqs. (15 1, (16), as one dimensional integrals. For example, Eq. (15 1 becomes 



Z±.d=7T~ <Pd,Bn-V<Pd,B<il , (17) 
27rw J C G 

where the integral on the right is a line integral over the curve that encloses the area of So, 
which we denote by Cq , and n is a unit vector normal to this curve in the direction of the 
axis (an identical integral over Cb was dropped, since <pd,B = on Cb)- Similarly, 



Zj_, q =—[ <p m ,Bn ■ ^<Pq,B dl . (18) 
' ic G 



TTU! 



The examples dealt with in this section are: (II) a small, flat iris (or horizontal slot) 
in a beam pipe, (lib) the same but with the design orbit shifted vertically, (12) a small 
rectangular iris, (13) a small elliptical iris, (14) a flat iris (not necessarily small) in a flat 
beam pipe. In all cases except lib, the design orbit is on a horizontal symmetry plane. A 
cross-section view of the geometries is sketched in Fig. [6] Dimension labels and the design 
orbit location are also shown. 

In the following subsections (except for subsection 14) we will assume that the pipe radius 
is very large compared to the iris aperture. We can use the free-space Green function for a 
particle vertically offset by z/o : 

G(x,y,y ) = -Hx 2 + (y-y ) 2 } . (19) 

From this Green function we obtain the potentials: 



2y 



Then using these potentials, we perform the calculations of Eqs. ujh, (U8 
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1 1 : narrow, flat iris 



lib: with shifted design orbit 
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12: small, rectangular iris 13: small, elliptical iris 





14: flat iris in flat beam pipe; longitudinal view 
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FIG. 6: Cross-section views of the iris/short collimator geometries considered in this section (for 
case 14 a longitudinal view is also given). Dimension labels are given, and the design orbit location 
is indicated by the red dot. In cases 11-13 the beam pipe aperture is specified as large compared to 
the iris aperture; in case 14 the beam pipe is flat, with vertical aperture 2b. For case 14 two colors 
are used for the boundaries as an aid in visualization. 



The error in transverse impedance introduced by the approximation of the free-space 
Green function is relatively small. For example, in the case of a round (i.e. cylindrically 
symmetric) iris of radius g within a round beam pipe of radius b it is known that the 
correction to the first order term is (g/b) 4 [S]. 



II. Flat Iris or Horizontal Slot in Large Pipe 



Consider a large beam pipe containing an iris with a horizontal slot of vertical size 2g, 
with g very small compared to the beam pipe size. We take the reference trajectory to be 
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the symmetry line, so the monopole term Zj_ iTO = 0. The dipole term is given by 

Z±,d= — — / tpd,B 0, 9) dtp J" B (x, g) dx , (21) 
nuj J dy 

and similarly for the quadrupole term. With the free-space potentials the integrals can be 
performed analytically, and after some algebra we obtain 

1 2 

Z±,d = Z± „ = — = , Z± = — ^ . (22) 
Log 1 tog 2 - 

We see that the dipole and quad terms are equal, and that the total result is the same as 
the leading order impedance for a round, thin iris of radius g [8]. Note that for a flat iris 
that is oriented vertically, we have Z± t d = —Z± j9 = \/{Log 2 ) and Z^ = 0. 



lib. Case of Shifted Design Orbit 

We can also find the impedance for the case the design orbit is shifted vertically by Ay. 
In this case the (transverse) monopole term Z± >m dominates the transverse impedance, and 
we can neglect the effect of Z± t d and Z± <q [see Eq. Q]. The transverse monopole term is 
obtained using Eqs. |5]), (|6j. Then the calculation proceeds in the same manner as before: 



Z 



_L,m 



1 



Vd,Fs{x,g- Ay) 



dtp* 



(x,g- Ay) 



+ VdM x i -9 ~ Ay) d<P ™> B ( x , -g - Ay) 



dy 



dx 



We find that the transverse impedance is given by 



(23) 



1 



1 



9 - Ay 



Ay 



(24) 



For small Ay, Z±_ m — 2 Ay /(tog 2 ), which is consistent with our earlier results. 



12. Rectangular Iris 



For the case of a rectangular iris with a small aperture of 2w by 2g (horizontal by vertical) 
the calculation follows the same procedure as before; however, the integration now is a line 
integral along the rectangular aperture. We perform integrals like 



J J-,d 



2 

TTLO 

The final solution is 

Z±,d = 

Z±,q — 

z±_ = 



<Pd,B( X >9) 



d(fid,j 
dy 



(x, g) dx 



<Pd,B(w,y) 



dtfd,j 
dx 



(w,y)dy 



(25) 



2 a + arccot(a) + a 2 arctan(a) 
TTLog 2 a 2 

2 a(—l + a 2 ) + (1 + a 2 )[— arccot(a) + a 2 arctan(a)] 
TTCog 2 a 2 {\ + a 2 ) 

4 a + (1 + a 2 ) arctan(a) 



tt Log 



l + a 2 



(26) 
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FIG. 7: For a rectangular iris in a beam pipe, the transverse impedances Zj_, Zx,q, normalized 

to 2/(ujg 2 ), as functions of (3 = (w — g)/(w + g). Plotting symbols give ECHO numerical results. 

where a = w/g. 

The results, when normalized to 2/uig 2 , are plotted as functions of /3 = (w — g) / (w + g) 
in Fig. [7] Note that when (3 = 1 the results agree with the results for the iris with the 
infinitely wide horizontal slot, given above. For j3 = — 1 (an infinitely high vertical slot) 
Z± lQ — —Zi_,d — —1/(luw 2 ) and Z± = 0. For the special case of a square aperture Z± = 
Z± lC i — 2(- + \)/g 2 and Z± tQ = 0. Note that the horizontal impedance is obtained from 
Eqs. ( |26[ ) by exchanging w and g. 

In Fig. [7] ECHO numerical results for Z± y d and Z±^ q are shown by the plotting symbols. 
We see good agreement with the results of our method. 



13. Elliptical Iris 

The impedance calculation for an elliptical iris with a small aperture, with axes w by 
g (horizontal by vertical), follows in a similar manner to the rectangular iris case. The 



equations that need to be solved are Eqs. (17), (18 1. We see that the elliptical case is more 
complicated in that it requires the formation of n • Vy 9i B and of the length metric on the 
ellipse of the iris, though the final solution is quite simple. We find that 

Z ±td = ~(l + ^), Z ± , q =~(l-^) , Z^=— 2 . (27) 

We find that for an elliptical iris Z\_ is independent of the size of the horizontal axis w\ 
As was true for the rectangular iris, the horizontal impedance of the elliptical iris can be 
obtained by exchanging w and g in the solution equations. In Fig. [8] we plot the analytical 
results and compare with ECHO numerical results. We see good agreement. 
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FIG. 8: For an elliptical iris in a beam pipe, the transverse impedances Zj_, Z± t 4, Z± tq , normalized 
to 2/(ujg 2 ), as functions of [3 = (w — g) / (w + g) . Plotting symbols give ECHO numerical results. 



14. Flat Iris in Flat Beam Pipe 

Consider a flat iris of height 2g centered within a flat beam pipe of height 2b. We begin 
this problem with the Green function between two parallel plates of aperture 2b |13j 



G(x,y,y ) = In 



cosh ff 



COS 



T~(y+ya) 
2b 



cosh ff 



lb 



(28) 



The calculation procedure is the same as for the free-space-to-flat-iris example, Example II. 
The potentials are 



The final solution is 
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In 



cosh ff 



Try 
2b 



cosh |f 



COS 



21, 



7r sin ■ 

2b cosh 2 §f - 

8b 2 



9 ttv 

2 + cosh + cos ^ ) cosh cos Sf 



(cosh 



2 7nc 
2b 



Z \ j — 



^ csc 2 (7ra)[27r(l — a) + sin(27ra)] 

2u>g 2 



J-,8 



Z\ = 



ira 



^9 



csc(7ra) [1 + 7r(l — a) cot(Tra)] 

a) + sin(7ra)] 



— csc 2 (W2)[7r(l 



(29) 



(30) 



where a = g/b. These curves are plotted in Fig. [9] The round case, with g and b, represent- 
ing, respectively, the radii of the iris and of the beam pipe, (Z±) roun d = 2(1 / g 2 — g 2 /b 4 )/uj [8], 
is also shown (the dashes). We note that Z± is always close to and larger for the flat than 
for the round case. 
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FIG. 9: For a flat iris with aperture 2g in a flat beam pipe of aperture 2b, the transverse impedances 
Z±, Zj_ : d, Z± :q as functions of a — g/b. Plotting symbols give ECHO numerical results. The round 
case, (Zx)round = 2(l/p 2 — g 2 /b 4 )/uj, is also shown (dashes). 

In Fig.[9]ECHO numerical results for Z±^ and Z± q are again shown by plotting symbols. 
We see basically good agreement with the results of our method. We should point out one 
subtlety however. For g/b = 0.8 the numerically obtained wake begins to take on some 
inductive character and the result for impedance begins to deviate from our analytical 
solution (particularly Z± <j)- We attribute this to the fact that, for the parameters used in 
the simulations, (6 — g)/b = a z ; as (6 — g)/b < <r z we are beginning to leave the optical 
regime. 



Longitudinal Impedance 

If we know the geometry of both the beam pipe and the iris we can also obtain the 



longitudinal impedance, Z\\ i ong . For the longitudinal impedance we begin with Eq. (Ill, 
which we convert to a line integral, and obtain 

2 r° dcp m , B f v . 

Z\\,long = / Vm.B—R [X, g) dX 

ncj dy 

sin 7rcv cosh irx , ( cosh ixx — cos 7ra \ , . 

In — dx , (31) 



cos 2na — cosh 2nx \ cosh ttx + cos 7ra 
th a = g/b. The last integral we evaluate numerically. The result, when normalized to 



4/c, is shown in shown in Fig. 10 For comparison, the round case impedance (with g and 
b representing, respectively, the iris and beam pipe radius), is also shown. We see that the 
longitudinal impedance in the flat case is always less than in the round case. 
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Z||,Long/(4/c) 




FIG. 10: For a thin flat iris with aperture 2g in a flat beam pipe of aperture 2b, the longitudinal 
impedance Z\\j. ng as function of a — g/b. The round case, (Z«) ro und = 41n(fc/g)/c, is also shown 
(dashes). 



IV. STEP-IN TRANSITION 

For the dipole component of impedance of a step-in transition we replace, in Eq. ([8| , Sc 
by Sb, and use Green's first identity to obtain 



Z\\.d = 



1 

2ttc 



V Cb J Sb 

tfd,Bn-X7tp diA dl+ fd,B^ 2 fd,AdS 

Cb J Sb 



(32) 



The first and third integrals are zero because fd,B = on boundary Cb, the second and 
fourth integrals cancel because the Laplacian is the same independent of region. 

The same kind of analysis shows that Z\\ q = in a step-in transition, and that Z\\, m = 
also in non-symmetric step-in transitions. Since specifics of the geometry were not used in 
our derivation, we conclude that the total transverse impedance is zero for a (short) step- 
in transition of any geometry. Similarly the longitudinal impedance of any short step-in 
transition is also zero in the optical regime (something that has been known to be true for 
the special case of a round step-in transition 0], [5]). We conclude that the longitudinal 
and transverse impedance of any step-in transition, in the optical regime, is zero! 

The impedance of a long collimator is just the sum of the impedances of a step-in and a 
step-out transition. Since the impedance of a step-in transition is zero, the impedance of a 
long collimator is the same as that of a step-out transition alone, the impedance of which 
we will study in the following section. 



V. STEP-OUT TRANSITION 



An example step-out transition is sketched in Fig. [TTJ (the particles move in the +z 
direction). The examples dealt with in this section are: (Tl) a flat beam pipe that transitions 
to larger flat beam pipe, (T2) a rectangular pipe that transitions to a large pipe, and (T3) an 
elliptical pipe that transitions to a large pipe. Cross-section views of the geometries (and 
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a longitudinal view in case (Tl) are sketched in Fig. 12 Dimension labels and the design 
orbit location are also shown. 




FIG. 11: A step-out transition. The particles move in the +z direction. 



For the impedance of a step-out transition we replace, in Eq. ([8]), Sq by Sa- The 
impedance function (in the dipole case) Zn d of a step-out transition can be written, be- 
ginning with Eq. as 



1 

2ttc 
1 

27TC 



/ iy^PdM) 2 dS-j X7ip diA -X7ip dtB dS 

J Sb J Sa 

^ (V<p d , B ) 2 dS- f (W VdA ) 2 dS 

Sb J Sa 



^Vd,A ■ ^{Vd,B - Vd.A) dS 



Sa 



By using Green's identity one can easily see that the last integral is zero. Thus 



Z 



1 

2nc 



I {V<p d , B ) 2 dS- [ {V^ dA ) 2 dS 

J Sb J Sa 



(33) 



(34) 



This result can be interpreted as 2/c times the static field energy of a dipole in region B 
minus that of a dipole in region A. This principle that the longitudinal impedance (here 
longitudinal dipole impedance) of a step-out transition in the optical regime is given by 
twice the difference in the static field energy in the two regions — was first elaborated for 
the round step-out transition (longitudinal case) by Heifets and Kheifets (Hj. The principle 
was then used for obtaining higher order wakes (azimuthal mode number m > 0) of (long) 
round collimators [7] , and also for specific 3D structures [5] . 

We see that a main difference in the calculation of the impedance of a step-out transition 
and of an iris is that, in the former case the surface integrals are performed over a region 
that includes the source charges of the potentials, and in the latter case the region of the 
charges is excluded. In the step-out transition case, however, it turns out that even though 
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Tl: flat step-out transition 
to flat beam pipe; 



t 

2 b 

i 



T 

2g 

i 



longitudinal view 



- • - - 2b 



T 

2 b 



T2: rectangular step-out 
transition 



T3: elliptical step-out 
transition 





FIG. 12: Cross-section views of the step-out transition geometries considered in this section (for 
case Tl a longitudinal view is also given). Dimension labels are given, and the design orbit location 
is indicated by the red dot. In cases T2-T3 the outgoing beam pipe aperture is specified as large 
compared to the incoming beam pipe aperture; in case Tl the outgoing beam pipe is flat, with 
vertical aperture 26. For case Tl two colors are used for the boundaries as an aid in visualization. 



the potentials diverge at the source charges, for the impedance two terms with the same 
divergence are subtracted, and the impedance is finite. 

By regrouping integrals and using Green's identity we can write the impedance of a step- 
out transition in terms of a simple line integral. The impedance function (for example, in 
the dipole case) can be written, beginning with Eq. d8j), as 



1 

2ttc 
1 

2ttc 



/ (V^d^s) 2 dS - V<fd,B ■V(<fd,A- (fid,B)dS 

J Sb — Sa J Sa 



tfd,Bn ■ Vfd.B dl- l <y5d,sV tp d ,B dS 

Cb—Ca J Sb — Sa 

I VdMV- • V(¥d,A ~ Pd,B) dl+ (Pd,B^ 2 (Pd,A - <Pd,B) dS 
J Ca J S a 



i 

2ttc 



tan' \7<f d ,Adl 



(35) 



To go from the second to the third form of the impedance in this equation, we use the fact 
that ifid.B = on Cb, that V 2 (pd,B = in the region Sb ~ Sa, and that ~S/ 2 {tfd,A — fd,B) = 
everywhere. 

However, we can obtain formulas for the impedance of a step-out transition that are even 



1G 



simpler. The dipole part of the impedance, for example, can be written as 



1 

2?rc 



C B 



Sb 



(fid,An ■ V(/3d,B dl 



C A 



: ij4 V ipd.sdS 



Sa 



(36) 



The line integrals arc zero again. Combining the remaining terms and substituting for V 2 y? 
from Eq. ([7]), we obtain 



2 d 

l±.d= lim — — [ip d .B(x, y) - (fid ,a(x,v)} 
x,3/-»0 u> ay 



(37) 



Both dipole potential terms diverge in the limit x,y —> 0; however, the divergences are the 
same, and Z± t d is finite. 

For the quad component of impedance of a step-out transition we replace, in Eq. (|9]), Sq 
by Sa, and use Green's first identity to obtain 



1 

2ttc 



fq,B^ 2 (Prn,BdS - 



Sa 



Vm,A^ 2 Vq,B dS 



(38) 



whe re the line integral terms are again zero. Substituting for the Laplacians from Eqs. 0), 
(101, using the relation \\m x ^ y ^[d 2 Lp„ h A{x, y)/dy 2 — 2ip qt A(x, y)] — 0, and going to we 
obtain 



Z ± , q = lim - [if q , B (x, y) - Vq.A{x, y)) 

x,y^0 UJ 



Similarly, the transverse monopole impedance is 

2 



Z±, m = lim - [(p d .B(x, y) - (p d .A(x, y)] 



(39) 



(40) 



the longitudinal impedance is given by 



,long 



lim - [(p m . B (x, y) - if m .A(x,y)} 

x,y—>0 C 



(41) 



We see that, once the Green functions are known, obtaining the impedances of step-out 
transitions (as well as long collimators) is a relatively straightforward matter. Whereas an 
equation of the form Eq. ( 34 1 can be interpreted to mean that the longitudinal impedance 
in the optical regime is 2/c times the static field energy of a charge distribution in region B 
minus that of one in region A, these simpler equations can be interpreted to mean that the 
longitudinal impedance is 2/c times the potential difference at the charge distribution in 
region B minus that at the charge in region A. 

As was the case with the irises/short collimators, for a transition into a large beam pipe, 
one also doesn't need to know details of the large beam pipe in order to calculate the leading 
order behavior of a transverse impedance. In the equations for impedance one can take the 
potentials in region B to be derived from the free-space Green function, given in Eq. (19 1. 
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Tl. Flat Transition 



In the case of a flat, symmetric, step-out transition going from aperture 2g to 2b, we 
substitute the potentials Eqs. (29 1 into Eqs. (37 1, (39 1. We obtain 



Z\ = 



2lo 



1 



1 

b 2 



(42) 



with Z± >q = 2-Z_L,d- We see thus that the transverse impedance of a flat step-out transition 
(or of a long, flat collimator) is a factor 7r 2 /8 times the transverse impedance of a long, 
round collimator, if we take the half-heights in the former case to be equal to the radii 
in the latter [BJ. In Fig. 13 we plot the theoretical dependence and compare with ECHO 
numerical results (the plotting symbols). We see that the agreement is very good. 



Z_L/(3i2/2t»]g2) 




FIG. 13: For a step-out transition from a flat pipe of aperture 2g into a flat pipe of aperture 26, 
the transverse impedances Z±, Z± t a, Zx,q as functions of a = g/b. Plotting symbols give ECHO 
numerical results for comparison. 



If we perform the longitudinal impedance calculation for the flat step-out transition we 
find that Z\\ i ong — 41n(6/g)/c, which is the same as for the round case, if we take the 
half-heights in the flat case to be equal to the radii in the round one. 



For a design orbit that is shifted by Ay from the symmetry axis, we find, using Eq. (40 1 



1 7rAy 1 7rAy 

- tan tan 

9 2g b 2b 



(43) 



Note that for a beam close to the wall at e.g. y = g, Z± j7n = 2/[ui(g — Ay)], which is twice 



as large is we found for the close-to-the-wall impedance in the case of an iris (see Eq. ( 24 ) ) 



T2. Rectangular Transition 



Consider a rectangular pipe of width 2w and height 2g that transitions into a large beam 
pipe, and a design orbit that follows the symmetry line of the rectangular pipe. The Green's 
function for Poisson's equation in a rectangular pipe has been obtained by Gluckstern, et 
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al [H]. Their result is: 1 



G(x,y,y ) 



n=l 



9 cosn 2 9 . nir nn 
ncosh^ Sm 2^ (y + ff) Sm 2^ (yo + - 9) 



2y 



ln|[x 2 + (2/-y ) 2 ] 
In [x 2 + (y - y ) 2 ] . 



sinh 2 25 + COS 2 ^(j/ + i/o) 

sin 2 £(y-2/o) 



4g 

sinh 2 f2 



(44) 



The last term contains the singularity. The other terms are everywhere finite, and the sum 
in the first term converges well. 

For our calculation we take region B to have the free space Green function, i.e. the same 



as the last term in Eq. (44 1. Thus we have no singularities left in our impedances, since 



they all involve the difference of potentials in the two regions. Our final result is 



3ujg 2 

7T 2 

6ujg 2 



1+24 E r 



c27rma 



m—1 

oo 



1-24V 2m ;\, 

L j ^ _|_ g7r(2m— l)a 

m—1 



(45) 



with a — w/g. We perform the sums numerically. 
■K 2 /(2ujg 2 ), are plotted as functions of /3 = (w — g)/{w ■ 



The results, when normalized to 
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(Z± = z ± . 



g) in Fig, _ 

Note that when (3—1 the results agree with the leading order behavior for the flat transi 
tion, given above. For (3 = — 1 (a step-out transition from an infinitely high vertical beam 
pipe) Zj_^ q = —Zj__ d = —n 2 / {(Slow 2 ) and Z±_ = 0. For the special case of a square beam pipe 
Z± = Z±,d = (0.697)7r 2 /(2uig 2 ), which is 86% of the result for a round pipe with radius g, 
A/(ojg 2 ); and Z± tQ = 0. The plotting symbols in the figure give ECHO numerical results, 
and we see good agreement. Finally, note that the horizontal impedance is obtained from 
Eqs. (451 by exchanging w and g. 



1 Note that there is a typo in their Eq. (5.11). 
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FIG. 14: For a rectangular step-out transition to a large beam pipe, the leading order dependence 
of the transverse impedances Zx, Z±,d, Z±, q , normalized to tt 2 / (2uig 2 ), as functions of f3 = (w — 
g)/(w + g). The plotting symbols give ECHO numerical results for comparison. 



T3. Elliptical Transition 



Consider an elliptical pipe of horizontal axis w and vertical axis g that transitions into a 
large beam pipe, and a design orbit that follows the symmetry line of the elliptical pipe. The 
Green's function for Poisson's equation in an elliptical pipe has been obtained by Gluckstern, 
et al [T3]. For the case w > g it is given by 2 



G(x,y,y ) 



- 4 E — 

71=1 

In [x 2 + (y - y f 



ReT„(^)ReT„(^f) ImT„(^) ImT„(^) 



cosh 711*0 



sinh uuq 



(46) 



with T n the Chebyshev polynomials of the first kind, d 2 = w 2 — g 2 , and uq — arccoth (w/g). 
[Note that G, as given here, is not zero but constant on the elliptical boundary, a fact, 
however, that does not affect our results.] The calculation procedure is the same as for the 
rectangular step-out transition. In the elliptical case we find that 



0d(x,y) 



(-irim [T 2m ^{^)} 



(x,y)=-4j2 



e 2(2m-l)u _ I 

(-l)™Re[T 2m (^)] 



O 4m-tio 



i 



(47) 



with signifying the part of the potential that does not contain the singularity. Note that 
the horizontal (x) impedance of a step-out transition is equal to the (vertical) impedance 
of the transition after it has been rotated by 90°. Thus to obtain the impedances for an 



Note that there is a typo in their version of this equation, their Eq. (4.36). 
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elliptical step-out transition with w < g we calculate the x impedances for the rotated case, 
following the analogous procedure to what we use for finding the y impedances. 
We finally obtain 



Z± : d — 



Z±.q — 



1G 



2m- 1 



Wfl 2(0> - 1) ^ ^ 



2m— 1 



32 



^ 2 (« 2 "i)S(g±i) 2m +1 ' 



(48) 



with a — w/g (valid for all a). We have simplified the results by using the relation e 2u ° — 
(a + l)/\a — 1|. The sums are performed numerically. The results are plotted as functions 
of j3 = (w — g)/(w + g) in Fig. |15| The plotting symbols in the figure give ECHO numerical 
results, and we see good agreement. We see that we obtain the expected results for the 
round case (f3 — 0), Z±^ = 4/(w<? 2 ) and Z±, q = 0; for the flat case (/3 — 1), Z± y d — 2Z±, q = 
7r 2 /(3cug 2 ). For the case of a horizontally inflnitesimally narrow elliptical pipe (/3 = —1) we 
find that Z±,d ~ ~Z± >q = 1/(luw 2 ), and that the total impedance Z± — 2/(ujg 2 ). Note that 
in this limit Zj_ has the same value as we obtained for any small elliptical iris in a beam pipe 
(see Eq. (27 1). Details of how all three limits are arrived at can be found in the Appendix. 
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FIG. 15: For an elliptical step-out transition to a large beam pipe, the leading order dependence 



of the transverse impedances Z±_ , Z± 



normalized to 7T / (2ujg ) , as functions of j3 = (i 



g)/(w + g). The plotting symbols give ECHO numerical results for comparison. 



VI. MORE COMPLICATED TRANSITIONS 



In this section we give two examples that are neither an iris/short collimator in a beam 
pipe, nor a step-in or step-out transition. The examples are: (Ul)-misaligned flat beam 
pipes, and (U2)-LCLS-type rectangular-to-round transitions. See Fig. 
perspective view of a pair of the LCLS transitions is also given in Fig. [T7] 



16 A cut-away 
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Ul : misaligned flat beam pipes 




U2: LCLS type rectangular-to-round transition 




FIG. 16: Cross-section view (from upstream end; left figures) and longitudinal view (right figures) of 
more complicated transitions: Ul-misaligned flat pipes, and U2-LCLS-type rectangular-to-round 
transitions. Dimension labels are given; the design orbit location is indicated by the red dot. Two 
colors are used as an aid in visualization. 




L 



FIG. 17: A matching pair of LCLS rectangular-to-round transitions in cut-away perspective view. 
The distance L between transitions is much larger than the catch-up distance for the nominal bunch 
length. 

Ul. Misaligned Flat Beam Pipes 

Consider first two flat beam pipes with thick walls and aperture 2g that are perfectly 
aligned and joined at z = 0. The design orbit lies in the horizontal symmetry plane. Now 
imagine shifting the z > pipe vertically by Ay (|Ay| < g) and the z < pipe by —Ay, 
and keeping the design orbit unchanged. Note that the resulting transition no longer has a 
horizontal symmetry plane. 

Let us sketch out the calculation of the transverse impedance Z± >m for this structure (for 
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Ay > 0). The potentials for this problem are given by 

dG 



<Pm,A{x, y) = G(x, y + Ay, Ay) 
<Pm,B{%, y) = G{x, y - Ay, -Ay) , cp d>B (x, y) 



¥>d,A(x,y) = ^{x,y + Ay,y ) 
oyo 



2/0= Ay 



dG 
dyo 



{x,y- Ay,y ) 



(49) 



2/o = -Ay 



with the flat pipe Green function G(x,y,yo) given by Eq. (28 1, but with b replaced by g. 
Note that, for this geometry, the aperture Sq is the intersection of Sa and Sb- 
Beginning with Eq. ^ , and using Green's identity, we obtain 



1 

2ttc 



(fid,Bn ■ ^<Pm,B dl 



C B 



:,bV tp m B dS 



Sb 



I ip d:B n-\7ip m:A dl+ / Lp d ,B^ 2 f m ,AdS 

JCc J So 



(50) 



The first integral above is zero, because tpd,B is zero on boundary Cb, and the second 
and fourth integrals cancel. We are left with the third integral, which implies a transverse 
impedance (see Eqs.[5j: 

/ ip d ,Bn ■ Vip m ,A dl 



2ttuj 



~- [ <Pd.B(x, 9 - Ay) dl ' f '™> A ( g] g - Ay) dx 
2-kw J dy 



(51) 



Note that the contribution of the integral at y = —g + Ay is zero because this is a boundary 
of region B. We obtain, finally, the analytical result (valid for either sign of Ay) 



Zx,m — 



1 

tog 



sgn(Ay) -7r l 



9 J 



cot 



7rAy 



7T esc ■ 



irAy 



[|Ay| < g] , (52) 



9 9 

with sgn(x) meaning the sign of x. We note that Z±^ m is positive for positive Ay, and that 
it is odd in Ay. Also, note that as Ay — > g, Z± >m — > 3/[ui(g — Ay)]. The result of Eq. 52 



(for Ay > 0) is plotted in Fig. 18 Plotting symbols give ECHO results; we sec excellent 
agreement with our results. The function, 



(Z±) 



approx 




(53) 



gives a good approximation to the impedance (see the dashes in the figure). 

Obtaining the longitudinal impedance Z\\ i ong for the misaligned pipe, one follows a similar 
procedure. In this case we find that the solution is given by 



'||, long 



1 

7TC 

2 cos a 

7T 



<fm,B(x,g- Ay) dif ™ ,A (x,g- Ay) dx 



In 



dy 

cosh x — sin a 
cosh x + sin 3a 



[Ay > 0] 



dx 

cosh x — sin a ' 



(54) 



with a — \~KAyjg. We solve the integral numerically. The result, for Ay > 0, is plotted 
in Fig. 19 Note that Z^ i ong is even with respect to Ay. The plotting symbols give ECHO 
results, and we see reasonably good agreement with our results. As a scale comparison the 
impedance of a round, step-out transition, from radius g to radius g + Ay, (Z^i ong ) roun d = 
41n(l + Ay/g)/c, is also given in the plot (the dashes). 
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FIG. 18: For the misaligned, flat pipe of aperture 2g, the transverse impedance Z± tm as func- 
tion of misalignment parameter Ay/g. Plotting symbols give ECHO results. The approximation, 
{Z±, m )a PP rox, is also shown (dashes). 




Ay/g 



FIG. 19: For the misaligned, flat pipe of aperture 2g, the longitudinal impedance Zui ong as function 
of misalignment parameter Ay/g. Plotting symbols give ECHO results. As a scale comparison the 
impedance of a round, step-out transition, from radius g to radius g + Ay, (Z\\) roun d = 41n(l + 
Ay/g)/c, is also given (dashes). 



U2. LCLS Rectangular-to-Round Transition 



In the LCLS undulator region there are 33 pairs of rectangular-to-round transitions. The 
rectangular aperture has horizontal width 2w — 10 mm by vertical height 2g = 5 mm; the 
round aperture has radius a — 4 mm. The axes of the two pipes are aligned. The transitions 
are abrupt. The bunch length in this region is 20 /im (rms). Thus our optical regime formulas 
are applicable. Note that these transitions are neither step-in nor step-out transitions. In 
the LCLS undulator region the longitudinal impedance is the more important one, and this 
is the one we calculate here. A numerical ECHO calculation has found that the longitudinal 
impedance of a pair of transitions (one rectangular-to-round and one round-to-rectangular 
transition) equals 1.21/c [9]. 
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Beginning with Eq. (11 1, and using Green's first identity, we obtain 



Z \\.long = -7T— / ^m.-B 11 ' V^ m ,A d£ . (55) 



Here Sq is the intersection of Sa and Sb- The integration path Cq follows the rectangular 
boundary at the top and bottom, and the circular boundary on the right and left (see 



Fig. 16 1. In the rectangular-to-round case [or we can say the rectangular-to-circular (rtc) 
case] with impedance Z\\ r t c region A is the rectangular pipe, region B is the circular pipe. 
In the circular-to-rectangular (ctr) case with impedance Z\\ ctr the regions are reversed. 
The circular monopole potential is given by ip m ,c — — ln[(£ 2 +y 2 )/afl; an d the rectangular 



monopole potential, by ip m ,r = G(x,y,0), with G given by Eq. (44). Then we have as 
rectangular-to-circular impedance 



z \\,rtc = — — / <Pm,c( x i9) ' (x,g)dx . (56) 



TTC J dy 

The contribution from the circular part of the boundary is zero, since ip m ,c is zero on this 
boundary. In the circular-to-rectangular case 



2a /• arctan (9/V a2 ~9 2 ) d(p m c 

Z\\,ctr — / (fim^r (a cos 9, a sin 9) — (a cos 9, a sin 9) d9 , (57) 



TTC J or 

where dtp miC /dr(x,y) = —2j\Jx 2 + y 2 . The contribution from the rectangular part of the 
boundary is zero, since ip m ,r is zero on this boundary. These integrals can easily be solved 
numerically, and the sums coming from ip m ^ r converge well. 

To do a small parameter study, let us keep the shape of the rectangular pipe fixed, with 
w = 2g, and let us vary a from g to w. In Fig. [20]we plot the results, giving the impedance 
of a rectangular-to-circular transition, Z^ rtc , of a circular-to-rectangular transition, Z^ ctri 
and the sum of one of each type, (Z\\) tota i, as functions of a/g. We see that each of the single 
transition curves goes to zero exactly where that transition becomes a step-in transition. 
For the actual design of the LCLS transition (a/g = 1.6), the rtc transition has 7.5 times 



the impedance of the ctr transition. In Fig. 20 the black dot gives (Z\\) to tai as obtained by 
ECHO, and we see good agreement: our calculation gives 1.24/c and the ECHO result is 
1.21/c. 



VII. CONCLUSIONS 



We have used a method, that we derived in a companion report [3], to find impedances in 
the optical regime, and applied it to various 3D beam pipe transitions that one encounters 
in vacuum chambers of accelerators. The method is applicable to high frequencies and 
transitions that are short compared to the catch-up distance. Our examples are of four 
types: an iris/short collimator in a beam pipe, a step-in transition, a step-out transition, 
and more complicated transitions. (Note that a long collimator with ends that are short 
transitions has an impedance that is the sum of the impedances of a step-in and a step- 
out transition, and is thus also included.) Most of our results are analytical, with a few 
given in terms of a simple one dimensional integral. We believe that all of our results are 
new. We have also compared (most of) our results with numerical simulations with the 
computer program ECHO, a finite-difference program that solves Maxwell's equations of an 
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FIG. 20: Longitudinal impedance for transitions of the LCLS rectangular-to-circular type, giving 
Z\\,rtc, Z\\,ctr, and their sum (Z^) to tai as functions of circular radius a. The rectangle width 2w — 4g. 
The ECHO result for (Zu)totah from Ref. [5], is given by the black dot. 

ultra-relativistic bunch within metallic boundaries of 3D geometry, and the agreement is 
excellent. Note that our method is a much simpler way of obtaining impedances than the 
simulations. 

We have focused on transverse impedances. For bi-symmetric (horizontal and vertical 
mirror symmetric) examples, a bunch moving at a small offset from the symmetry axis will 
excite a transverse impedance composed of both a dipole and a quadrupole component. For 
such problems we give both components. The iris/short-collimator-in-a-beam-pipe examples 
we solve include irises with small aperture that are (II) flat, (12) rectangular, (13) elliptical; 
also included is (14) a flat iris (not necessarily small) in a flat beam pipe. An interesting 
result is that the vertical impedance of an elliptical iris is independent of the horizontal axis 
of the ellipse. 

For a step-in transition of any shape we find that all impedances (transverse and lon- 
gitudinal) are zero. For a step-out transition (which in the optical regime has the same 
impedance as a long collimator) we give the solution for (Tl) a flat step-out transition to 
a flat beam pipe, (T2) a rectangular step-out transition, and (T3) an elliptical step-out 
transition. We find, for example, that the transverse impedance of a flat, long collimator in 
a flat beam pipe is 7r 2 /8 times the impedance of the inscribed round, long collimator in a 
round beam pipe. The more complicated transitions examples we solve are (Ul) misaligned 
flat pipes and (U2) the LCLS rectangular-to-round transitions. 

The method of Ref. [3] is powerful; it allows one to calculate the impedance in the optical 
regime of a truly large class of transitions. We have demonstrated this with a small number 
of relatively simple examples, compared to what is possible. 

Acknowledgements 

We thank CST GmbH for letting us use CST MICROWAVE STUDIO for the meshing for 
the ECHO simulations. We also thank B. Podobedov for running a GdfidL test simulation 
for us. 



2G 



Appendix: Limiting Values of the Impedance of Elliptical Transitions 



The solution for the elliptical step-out transition with axes w by g (horizontal by vertical) , 
Eqs. [48] can be written as 

4(1 -(3) 2 ^ 2m -I 



_ 8(1 -/3) 2 m 
Z± ' q cog 2 ? ^ 0-*™ + 1 ' (58) 

where f3 — (w — g)/(w + g). We derive here the limits for /3 = (a step-out transition from 
a round beam pipe), /3 = 1 (from an elliptical pipe with infinitesimal height), and f3 — — 1 
(from an elliptical pipe with infinitesimal width). 



The Limit for (3 = 

In this limit, only the m = 1 term contributes to Z± d and no term contributes to Z± q . 
We obtain the round step-out transition results 



4 



2 



Z Uq = . (59) 



The Limit for = 1 

Let us consider the dipole impedance first. For (3 — 1 — e, with e a small positive number, 
the sum peaks when 2me ~ In m, i e. for a large value of m. Thus the sums can be replaced 
by integrals, ^2 m — » / dm. Changing variables to x = /3~( 2m_1 ), dx = — 2(ln/3)a; dm, we 
obtain 

^-^lo^^hWA (^lyv (60) 

The integral equals 7r 2 /6, and 

«w = • (61) 

The calculation for Z±^ q is similar. Our final result is the same as for the flat pipe step-out 
transition 

2 2 2 

^ d= 3^' Z ^ = £^' Z±= 2^- (62) 
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The Limit for /3 = — 1 



Let 7 



-/?. Then our equations become 
4(1 + 7 ) 2 
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m 

-2m _ 



(63) 



To find the leading order behavior of Zj_ we need to go to second order in the calculation. 
To do this we will use the Euler-Maclaurin formula relating sums to integrals |15j 



/(I) 



f(x)dx--f'(l) , 



(64) 



where ' denotes taking the derivative of a function. The formula is valid if the sum converges. 
The approximation is good if \f"'(x) \ dx is small. 

Let us consider first the dipole part. We want the solution for 7 — 1 — e, with e a small, 
positive parameter that we, in the end, let go to zero. For the dipole part 



f(m) 



16 



2m- 1 



Ujg 2 7 -(2m-l) + I 



(65) 



There are three terms on the right side of Eq. [64] that we need to calculate. The first term 
/(l)/2 = A/{uig 2 ). The integral term, done like before, is 



f(x) dx = lim 



1 



7->l ujg 2 (hi7) 2 J 1 / 1 (x+ 1) x 



In x dx 



(66) 



Note that for this order of calculation, in the integral on the right, the lower limit of 
integration is 1 /7 instead of 1 . This integral 



00 In x dx tt 2 

r — = m 

/ 7 (x + 1) x 6 



1 + 7 



In 7- i(ln 7 ) 2 +Li 2 (-l/ 7 ) 



(67) 



with Li 2 (a;) the polylogarithmic function of order 2. We combine this result with the terms 
in front of the integral in Eq. [66| expand around 7 = 1, and substitute 7 = (g + w) / (g — w); 
we find that the second, integral contribution to Eq. 64 equals 
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The third term in Eq. 64 — /'(1)/12 = — 4/(3g 2 ). We sum all three contributions to 
obtain Zj_ ^. Exactly the same technique is used for Zj_ q (note, however, that the equation 



corresponding to Eq. 
obtain 



(>(> 



will have an integral with lower limit I/7 2 , not 1/7). We finally 
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(68) 
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The loading order behavior of these impedances, in this limit, is Z± :< j = —Z±. q = it 2 /(6ojw 2 ). 
Summing the two impedances together, we find that the total impedance equals Zj_ = 
2/{ujg 2 ). This value of total impedance is the same as we found for any small elliptical iris 
in a beam pipe. This result appears to agree with the numerical calculation of the original 



sums, which is plotted in Fig. 15 
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